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Abstract 

■ The presence of gas in dusty debris disks around main-sequence stars is 
| ^ ! reviewed. We present new observational results for the most prominent rep- 
resentative of the class, viz. the southern naked-eye star 8 Pictoris. The spa- 

■ tial and spectral distribution of observed atomic lines from the disk around 
the star is reproducable by a Keplerian rotation model to a high degree of 
accuracy. The expected velocity dispersion due to radiation pressure in reso- 

■ nance lines is not observed. Modeling the motion of different atomic species 
Q\ . under the influence of gravity, radiation pressure and gas friction leads to the 
^ | conclusion that an underlying decelerating component must be present in the 

disk. This braking agent is most likely hydrogen, with inferred average den- 
. sities nji > 10 6 cm -3 . This could support the observational result of Thi 

et al. (2001) which indicated the presence of appreciable amounts of H2 
i-£h ■ around the star 3 Pic. 

Oh; 

6 ■ 
a: 

ctf ■ 

>• ! One of the major highlights of the IRAS mission was the discovery of the 'Vega- 

I excess' phenomenon exhibited by a number of main-sequence stars (Aumann et 

al. 1984), i.e. emission in excess over the stellar photospheric value at long wave- 
lengths (> 10/im). Already from the very beginning it was recognised that this 
excess is re-emitted stellar radiation by a surrounding cloud of dust particles. The 
image of the southern Vega-excess star (5 Pic in the optical by Smith & Terrile 
(1984) confirmed the earlier notion that the spatial distribution of the dust is highly 
flattened, i.e. in a disk shaped configuration rather than in a spheroidal cloud. Re- 
cent reviews of j3 Pic and similar systems - the dusty debris disks - include those 
by Artymowicz (2000), Lagrange et al. (2000) and Zuckerman (2001). 

'Collaborators on this project include Pawel Artymowicz, Alexis Brandeker, Goran Olofsson 
(Stockholm Observatory, Sweden), Malcolm Fridlund (ESA-ESTEC, The Netherlands) & Taku 
Takeuchi (University of California, Santa Cruz, USA). However, the responsibility for this article 
rests entirely with the author. 



1 Introduction 
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1.1 Dusty Debris Disk Evolution 

Besides Vega (a Lyr) and j3 Pic, other prominent members of the class are Foma- 
lhaut (aPsA), eEri and HR4796A. The excess emission of these objects seems 
to outline an evolutionary trend (Spangler et al. 2001, Fig. 1), in the sense that 
the amount of emitting particulate matter in a given object declines with time. 
Yet younger objects would occupy the upper left corner of the diagram (with 
considerable spread due to increasing optical depth) and it seems reasonable to 
assume that the dusty debris disks (somehow) are linked to those of the pre-main- 
sequence (PMS) phase. PMS-disks, with ages ^ lOMyr, are essentially gaseous 
and are b elieved to provide the material necessary for the formation of planetary 
systems ( Calvet et al. 200ol) . Observed asymmetries and inhomogeneties in the 



dust emission from debris disks have been taken to provide evidence for dynam- 
ical signatures of disk-planet interaction, (e.g., Ozernoy et al. 2000, Augereau et 
al. 2001, Wyatt & Dent 2002, Kuchner & Holman 2003). 

1.2 Planets in Dusty Debris Disks 

Our knowledge regarding the incidence of known planets (gas giants) around de- 
bris disk stars is very limited. Heavy observational bias is at work here, though. 
Radial-velocity studies are biased toward relatively late spectral types, whereas 
the most prominent debris disk stars are of earlier spectral types. In addition, 
radial-velocity studies have focused on fairly old stars, whereas disks around stars 
much older th an about 400 Myr w ould generally not have been picked up in the 
IRAS survey (Habing et al. 1999h . Three planetary systems, discovered with the 



radial-velocity technique, have been claimed by Trilling et al. (2000) to possess 
dust disks. The observational evidence presented by Jayawardhana et al. (2002) 
(see also Schneider et al. 2001) argues strongly against any substantial amounts of 
dust around 55Cnc, whereas the other two systems, viz. pCrB and HD 210277, 
still need confirmation. Pantin et al. (2000) reported the detection of a circumstel- 
lar disk around i Hor. However, this result has not yet been published in a refereed 
journal. 

The circumstantial evidence for planet(s) around dusty debris stars is enor- 
mous. To give just a few examples, the case for Vega has been presented by 
Wilner et al. (2002), for Fomalhaut by Holland et al. (2003) and for (3 Pic by 
Wahhaj et al. (2003). For the chromospherically active star eEri, Hatzes et al. 
(2000) presented a rather noisy radial velocity curve (for more recent data, see 
http://exoplanets.org/esp/epseri/epseri.shtml) and the confirmation of the planet(s) 
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would need a longer baseline in time. 



1.3 Widespread Gas in Dusty Debris Disks 

If the evolutionary hypothesis is correct, then a similar apparent decrease in the 
gas content may be expected and one could hope to find objects being currently in 
a transitional phase between gas-rich PMS- and gas-poor debris disks. This could 
open up the possibility to directly witness the end phases of planet formation and 
the early stages of planetary system evolution. Appreciable amounts of gas could 
still be present in these systems. 

Alternatively, in older systems, where the dust debris is most likely of more 
recent origin, the study of the dust production mechanism(s) could benefit from 
complementing observations of the accom panying gaseous compo nents, as for 



example in the case of solar system comets dCremonese et al. 20021) . 

In either of these scenarios, the detection of sharp spectral features originating 
from the gas would permit the direct observational study of the disk kinematics 
and provide input to theoretical models of the debris disk dynamics. 



2 Searches for Molecular Gas 

Recalling that the discoveries were made at IRAS wavelengths, it was clear al- 
ready from the very beginning that the emitting dust grains had to be at relatively 
low temperatures (a few tens to some hundred Kelvin). The emission being opti- 
cally thin, these observations led dire ctly to estimates o f the total amount of dust 



(particle sizes < A) in these systems (Hildebrand 1983). Disregarding for a mo 



ment the many details which in reality enter these estim ations such as, e.g., t he 
correct description of the dust absorption coefficient(s) dBeckwith et al. 200oT) - 
the mass is inversely proportional to this quantity - inferred lower limits to the 
dust mass are some fraction to some tens of a lunar mass (M Mo on)- For compar- 
ison, the Minimum Mass Solar Nebula amounts to roughly 5 x 10 4 M Mo on> from 
which it is clear that, for any reasonable assumption about the local gas-to-dust 
mass ratio, any solar system type of planet formation must have been completed 
by the Vega excess stage. 

If the gas-to-dust mass ratio in these systems were comparable to that of the in- 
terstellar medium (ISM, ~ 160), detectable amounts of (primarily molecular) gas 
could be expected. In fact, carbon monoxide (CO) has frequently been observed 
through disk absorption against the ultraviolet stellar continuum of (3 Pic (Roberge 
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et al. 2000 and references therein). In an optical absorption study, Hobbs et al. 
(1985) put an upper limit on the amount of CH + in the (3 Pic disk. 

Emission line studies, mostly in the sub-/millimeter regime, have so far been 
essentially unsuccessful (Liseau 1999 and references therein). Molecules searched 
for in emission include CO, CS, SiO and HCO + from the gound and H 2 and O2 
(the SWAS Team) from space. In addition, the detection of rotati onal H? emissio n 
in the mid-infrared (ISO-SWS) from /5Pic has been reported (Th i et al. 20011) . 
with an estimate of the total mass of the order of 50 Af®. However, Lecavelier des 
Etangs et al. (2001) were unable to confirm these amounts of H 2 through absorp- 
tion measurements in the far ultraviolet (FUSE) and the matter seems at present 
inconclusive (but see below, Sect. 3.1.2). 

Various explanations have been offered to explain the observed low levels of 
molecular concentration in dusty debris disks, particularly what regards CO, be- 
ing ubiquitous in the ISM and in solar system comets. These include 'abnormally' 
low gas-to-dust mass ratios, 'abnormally' low gas phase CO abundances (rela- 
tive to H 2 ), photodissociation of CO and, simply, non-existence (early dissipation 
through stellar wind and/or consumption during planet building). None of these 
explanations is flawless and the resolution of this issue will have to await future 
observational and theoretical improvements. 



3 Searches for Atomic Gas 

Observational studies of atomic gas, either ionised or neutral, in dusty debris disks 
have traditionally focused on line absorption in a very limited number of objects. 
Although very sensitive to column densities, these studies are generally incapable 
of pin-pointing the regions of line absorption along the line of sight. In the fol- 
lowing, a number of observational studies are discussed, which were aiming at the 
detection of emission from atomic gas. 



3.1 Atomic Gas Emission from the (3 Pic Disk 

The by-far best studied object among the dusty debris disks is the (3 Pic system. 
The attempt to directly measure, throu gh 21cm line emissi on, the amount of 
atomic hydrogen gas was unsuccessful (IFreudling et al. 1995b . The presence of 
widespread ionised hy drogen in the disk is not expected and B aimer line emission 
has not been detected dBrandeker et al. 2002r) . Observations with the ISO-LWS 
may have detected feeble [C 11] 157 /im emission toward (3 Pic (IKamp et al. 20 03). 
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Few details about the circumstellar disk could be learned from these spatially and 
spectrally under-resolved observations, though. 

The resolution issue was greatly improved upon by using a large optical tele- 
scope, equipped with an Echelle spectrograph and a long slit, leadin g to the dis- 



covery of widespread sodium gas in the disk (lOlofsson et al. 20011) . These ob- 
servations provided finally conclusive evidence that the star (3 Pic is indeed sur- 
rounded by a disk in Kepleri an rotation, seen nearly edge-on (Fig. 2). 



Follow-up observations (Brande ker et al. 20031) covered the optical spectrum 



(0.3 to 1.1 /im) at greater resolving power and, in particular, covered a larger area 
with the spectrograph slit positioned both along and perpendicularly to the disk- 
midplane (Figs. 2 and 3). Projected slit widths were 6 AU and slit lengths 154 AU, 
covering the disk out to 300 AU in both the northeast (NE) and southwest (SW) 
parts of the disk. Perpendicular slits were positioned at approximately 60 AU 
and 120 AU on both sides from the star, extending to 80 AU on either side of the 
midplane (Fig. 3). 

Nearly 80 emission lines from the disk were detected, all from non-volatile 
elements including Call, CrI+II, Fel, Nal, Nil+II, Til+II. In addition, a dozen 
lines lack so far a clear identification. Neutral iron, with more than 50 strong 
lines, dominates the disk spectrum. Measurements of the line emission along the 
perpendicular slits revealed that the gas disk is not perfectly plane (Fig. 4). Inside 
±3", the disk appears tilted at an angle of about 5°, very much the same as inferred 
for the dust disk by Heap et al. (2000). Further out, the gas disk is flaring at an 
opening angle of about 10° (H/r ~ 0.2). 

The two disk halves display significant asymmetries and small scale structure 
in their line emission, as illustrated by NaID 2 in Fig. 5. The NE disk extends 
much farther out, in fact to the limit of our observations (cf. Fig. 2). In contrast, 
further in, the SW side is brig hter. These featur es can also be discerned in the dust 
disk, albeit at lower contrast ( Heap et al. 2000l) . 



3.1.1 The Kinematics of the Atomic Gas Disk 

As a complete surprise came t he observed small ve locity dispersion of the line 



emission from the atomic disk (lOlofsson et al. 20011) . In the absence of a major 



component of volatile gases, an excess of some hundred of km s~ x over the local 
Keplerian value of the radial velocity would be expected due to radiation pressure 
on the neutral sodium atoms. As different species will be differently sensitive to 
radiation pressure, their spectral lines should exhibit different widths. Below, the 
radiation pressure coefficient is calculated for various atoms/ions observed in the 
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(5 Pic disk. 

For a single transition i, the radiation pressure coefficient expressing the 
relative importance of radiation pressure and gravity, can be written as 

' f iniogs* A + ' 

± grav ±u -^atom 

where the flux density Fa refers to the stellar photosphere at the wavelength 
A of the line and the gravity is given by the surface logg* and the atomic mass 
number A atom of the ion, expressed in amu, and where 

7rr e /i = -^A 2 4d • (2) 

r c is the classical electron radius, /; the absorption oscillator strength and the 
other symbols have their usual meaning. For a particular atom/ion, the radiation 
pressure coefficient is (3 = J2i A- Both F grav and F rad are inverse square laws so 
that their ratio is independent of the distance to the star. 

In Table[H /5-values for a number of atoms/ions are presented. These are all 
for transitions from the ground, as the disk gas is too cool to maintain significant 
excitation in the higher states. F or the photospheric spe ctrum the Next Generation 
Atmosphere Models were used ( Hauschildt et al. 19991) . which is shown in Fig. 6. 



Two values of (3 are given per resonance line. j3 min corresponds to the scattering of 
a stellar photon by an atom which is essentially at rest (with respect to the radial 
velocity component), i.e. the source is the core of the photospheric absorption line. 
As the atom picks up speed, a greater number of source photons becomes available 
from the wing of the photospheric line and /3 max refers thus to the neighbouring 
stellar continuum. In addition, as we did not spin up the model atmosphere to the 
rotation velocity of the star, these two j3 values respresent truly strict limits to the 
realistic case. Obviously, the differences can be quite dramatic, as for e.g. the H 
and K lines of Call (Table[T]). In contrast, the photospheric NaD absorption is 
quite shallow and, hence, /3 m i n ~ f3 m &x so that the sodium atoms should be sailing 
at their terminal speed at velocities much higher than what is actually observed 
(Olofsson et al. 2001, Brandeker et al. 2003). 



3.1.2 The Dynamics of the Disk Gas 

Also quite obvious from Tabled] is the fact that lines originating from different 
species and/or different ionisation stages should behave quite differently. For in- 
stance, the lines of singly ionised chromium should reflect only Keplerian rotation, 
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Table 1 : Radiation pressure coefficients a (3 for the j3 Pic gas disk; ions in alpha- 
betical order 



Ion Vacuum 

Wavelength (A) 



^4ui 
(s^ 1 ) 



# T 



1 

max 



Cal 

Sum (9) c 
Call 

Sum (4) 
CrI 

Sum (18) 
CrII 
Sum (5) 
Fel 

Sum (40) 
Fell 

Sum (14) 
HI 

Sum (10) 
Nal 

Sum (32) 
Nil 

Sum (36) 
NiE 
Sum (25) 



4227.918 

3934.777 
3969.591 

3579.705 
3594.507 

8002.280 

2484.021 
3861.006 

2382.765 
2600.173 

1215.6682 
1215.6736 

5891.583 
5897.558 

2290.690 
2320.747 

1751.910 



2 .18e+8 

1.47e+8 
1.4 e+8 

1.48e+8 
1.50e+8 

1.0 e-1 

4.9 e+8 

9.7 e+6 

3.8 e+8 
2.2 e+8 

6.265e+8 
6.265e+8 

6.22e+7 
6.18e+7 

2.1 e+8 

6.9 e+8 

4.8 e+7 



35 
35 
2.0 
1.0 
3.1 
3.8 
3.1 
14 

1. e-7 


0.16 
1.3 
3.7 

0.03 
0.4 
0.6 

3.e-5 

2. e-5 


160 
80 
240 
0.06 
0.06 
0.5 
0.06 




270 
270 

95 

34 
130 

16 

13 

69 

1. e-7 


11 
20 
32 

7 

5 

16 

2. e-2 
l.e-2 



170 
83 
254 
2.3 
12 
20 
0.2 
0.2 



Til 

Sum (40) 



3636.499 8.04e+7 
3982.887 3.76e+7 



2.0 
1.6 
14 



10 
10 

77 



Till 
Sum (7) 



3384.740 1.09e+8 
7 



3.0 
3.9 



12 
19 



a NIST atomic data; m odel is T cfr = 8000 K, \ogg = 4.5, logZ/Z G = 0.0 
( Hauschildtetal. 1999h . 

b (3 = Frad/^grav an d ™ and max values refer to the central line absorption and 

neighbouring continuum, respectively. 

c Number of around state transitions for total 8. 



whereas neut ral chromium should po ssess a significant radial component. This is 



not observed dBrandeker et al. 20031) . In fact, regardless of their j3 value, none of 
the species shows significant radial velocity excess over the Keplerian value. 

A possible explanation could be the existence of considerable amounts of, yet 
to be identified, 'quiescent' material, braking the radiatively accelerated gas to 
the observed relative velocities. An obvious candidate would be hydrogen, being 
hard to get at observationally. The star produces virtually no flux below the Lyman 
limit or at the Lyman transitions of H I (see Fig. 6). Any hydrogen gas would thus 
only follow the Keplerian rotation and be radially stationary and, as such, could 
act as the braking agent (see Table[l]). 

We apply a simplified analysis, including gravity, radiation pressure and gas 
'friction', to examine this possibility. The interaction of the gas with the dust 
particles in the disk is neglected (only big grains with /3 < 1 and strong coupling 
to the gas could decelarate the atoms). The equation of motion of the atom reads 



^ ,j -fgrav "I - F?aA Ffric • (3) 



dv 

dt 

Letting F ra d = j3 -F grav and assuming Ff ric — Cv, this reduces to 



m^ = (j3-l)F psr -Cv . (4) 
dt 

The radial solution to Eq. 4 is 

v = v oo + (v -v oo )e- ct / m (5) 
with v = v(t = t ) and, as t — > oo, 

^oo -^grav • (6) 

For collisions of neutrals with neutral hydrogen atoms, we approximate the 
friction coefficient (full momentum deposition) by 

Cnn = vr Oq m H n u v n = a 2 (8 vr k m H ) 5 n H , (7) 

where a is the Bohr radius. For ions colliding with neutral hydrogen atoms, 
Beust et al. (1989) proposed (hyperbolic orbit approximation) 

C in = 7tblm H n H v u = (it a H m u /e )^ qn R , (8) 
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where 6 is the impact parameter, a H is the polaris ability of hydrogen and q 
is the electric charge of the ion. The characteristic time scale to approach is 
given by the e-folding time 

T = m/C . (9) 

For Eq. 6 to be valid r needs to be shorter than the time scale for significant 
changes in the disk (e.g., |dF grav /F grav |~ 2v 00 r/r <C 1), a requirement which 
was verified a posteriori to be generally fulfilled. 

In Fig. 7, the asymptotic velocity Voo is shown as a function of the radial dis- 
tance from the star, r, for three atomic species having different atomic masses and 
very different values of [3 (Table[T]). The disk structure is approximated by power 
laws, e.g. n oc r p , where the density law was varied, but the same temperature 
distribution was used throughout [power law exponent —0.3 and normalisation at 
26 AU, i.e. T(26 AU) = 110 K]. Density distributions are shown for three exam- 
ple values of p, i.e. p = (constant density), p— —1.5 and p = —2.5 ('standard' 
density distribution). 

For the neutral species 23 Na I, an average value of (3 = 250 was used and in 
Fig. 7 the normalisation values at 26 AU of the hydrogen density [n(26AU) in 
cm" 3 ] are 10 6 for p = 0, 10 7 for —1.5 and 10 8 for —2.5. As an illuminating 
example, for the ion 40 Ca II, the minimum value /3 m i n = 3 was selected and the 
corresponding n(26 AU) are 10 4,5 forp = and 10 5 for both —1.5 and —2.5. As 
the radiation pressure coefficient of 52 CrII is merely a tiny 4 x 10" 7 , only the 
constant density case [n(26 AU) = 10 5 cm -3 ] is shown as an illustration. 

Whatever the actual gas density distribution in the (3 Pic disk might be, it is 
clear from these numerical experiments that average densities would need to be in 
excess of 10 6 cm" 3 in order to meet the limits set by the observations: the observed 
spatial distributio n of the line profiles a re reproduced to within 0.1 kms" 1 by a 



Keplerian model (Brand eker et al. 20 03). If this gas were roughly spherically dis- 



tributed, as hinted at by the observations perpendicular to the disk, the total mass 
would be in the neighbourhood of Mh ~ 0.02 r 3 00 /x 2 Mj upiter (rioo is in units 
of 100 AU; fi2 is the molecular weight in units of one particle per two hydrogen 
nuclei; n 6 is n H in units of 10 6 cm" 3 ). The implied column density of hydrogen 
Ah ~ 2.5 x 10 21 n 6 r 100 cm -2 would suggest this ga s to be molecular, as th e upper 



limit on atomic hydrogen is a few times 10 19 cm -2 dFreudling et al. 19951) . It fol- 
lows that such an H 2 component would be consistent with the observations of Thi 
et al. (2001) for r wo ~ 2. Confirmation of these results by re-newed observations, 
with e.g. SIRTF, are thus eagerly awaited. 
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Significant amounts of H 2 around (3 Pic are hard to reconcile with the absence 
of detectable H 2 absorption along the line of sight toward the star, at the same time 
as CO absorption has clearly been detected. If real, the existence of such relatively 
large amounts of gas in the (3 Pic disk would have far reaching implications for th e 
gas-dust dynamics and the evolution of the disk (Ta keuchi & Artvmowicz 20011) . 
However, as already pointed out, this gas component would be insufficient to build 
an analogue to the solar system and may simply represent 'left-overs' from a re- 
cent planet formation epoch. 



3.2 Search for Atomic gas in Other Systems 

To this end, we have also searched the disks around HR4796A and eEri for 
atomic gas. These searches have been unsuccessful. The reasons for this fail- 
ure are possibly quite different, though. If the amount of gas would scale with the 
amount of dust, the dustier disk around HR 4796A should exhibit an even clearer 
gas signature than the (3 Pic disk. However, HR 4796A is nearly four times more 
distant and the scatteri ng disk extends to less than one arcsec from the central star 
( Schneider et al. 1999b . presenting a most unfavourable contrast case for ground 
based observations. 

e Eri on the other hand is more than six times closer to the Earth than j3 Pic and 
its circumstellar disk subtends several tens of arcsec in the sky (at sub-millimeter 
wavelengths). It should thus, in principle, be easier to observe. The disk con- 
tains merely one percent of dust compared to the /3Pic disk (see also Fig. 1). 
Again, scaling the gas emission with the dust emission would imply two orders of 
magnitude lower intensities, which would have escaped detection at our present 
sensitivity. 



4 Conclusions 

Below, the main conclusions of this work are briefly summarised. 

• Searches for molecular gas in dusty disks around main-sequence stars have 
so far been unsuccessful, with the possible exception of the j3 Pic system, 
toward which Thi et al. (2001) reported the discovery of substantial amounts 
of H 2 gas. 

• Searches for atomic gas in a few of these systems h ave resulted thus far i n 



clear detections of line emission from the (3 Pic disk (Olofsso n et al. 20011) . 
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Follow-up observations revealed the presence of a large number of spectral 
lines from neutral and singly ionised non- volatile species (B randeker et al. 20031) . 



• The results from spatio-spectral fitting of these lines are consistent with gas 
orbiting the star at Keplerian velocities. Deviations from Keplerian motion 
are less than 0.5 km s _1 at the 5cr level. 

• The analysis of these motions invoking gravity, radiation pressure and gas 
friction on the atoms/ions indicates that some kind of decelerating material 
must be present in the disk. 

• Identifying this braking agent with hydrogen gas leads to the conclusion 
that average hydrogen densities are nn > 10 6 cm~ 3 and that the hydrogen 
is most likely molecular. 

• Without stretching parameters too far, this result could be in support of the 
result by Thi et al. (2001), i.e. that the present amount of gas around (3 Pic 
is of the order of 50 M®. 

Acknowledgements 

I wish to thank the organisers for inviting me to this very stimulating conference 
and Alexis Brandeker for Figures 2 through 5 and for valuable comments on the 
manuscript. 

I'd like also to mention that I am convinced that the collaboration between 
ESA and NASA will result in a highly successful Darwin/TPF mission. 

References 

[Artymowicz 2000] Artymowicz P., 2000, SSR, 92, 69 

[Augereau et al. 2001] Augereau J.C., Nelson R.P, Lagrange A.M., Papaloizou 
J.C.B. & Mouillet D., 2001, A&A, 370, 447 

[Aumannetal. 1984] Aumann, H.H., Beichman C.A., Gillett F.C., et al., 1984, 
ApJ, 278, L 23 



11 



[Beckwith et al. 2000] Beckwith S.V.W., Henning T. & Nakagawa Y., 2000, in: 
Mannings V., Boss A. P. Russell S.S. (eds.) Protostars and Planets IV, University 
of Arizona, p. 533 

[Beust et al. 1989] Beust H., Lagrange-Henri A.M., Vidal-Madjar A. & Ferlet R., 
1989, A&A, 223, 304 

[Brandeker et al. 2002] Brandeker A., Liseau R., Olofsson G. & Fridlund M., 
2002, in: Caroff L. & Backman D. (eds.), Debris disks and the formation of 
planets: a symposium in memory of Fred Gillet, ASP, in press 

[Brandeker et al. 2003] Brandeker A., et al., 2003, in preparation 

[Calvet et al. 2000] Calvet N., Hartmann L., & Strom S.E., 2000, in: Mannings 
V., Boss A. P. Russell S.S. (eds.) Protostars and Planets IV, University of Ari- 
zona, p. 377 

[Cremonese et al. 2002] Cremonese G., Huebner W.F, Rauer H. & Boice D.C., 
2002, AdSpR, 29, 1187 

[Freudling et al. 1995] Freudling W., Lagrange A.M., Vidal-Madjar A., Ferlet R. 
& Forveille T., 1995, A&A, 301, 231 

[Habing et al. 1999] Habing H.J., Dominik C, Jourdain de Muizon M., et al., 
1999, Nature, 401,456 

[Hatzes et al. 2000] Hatzes A.P, Cochran W.D., McArthur B., et al., 2000, ApJ, 
544, L 145 

[Hauschildt et al. 1999] Hauschildt PH., Allard F. & Baron E., 1999, ApJ, 512, 
377 

[Heap et al. 2000] Heap S.R., Lindler D.J., Lanz T.M., et al., 2000, ApJ, 539, 435 

[Hildebrand 1983] Hildebrand R.H., 1983, QJRAS, 24, 267 

[Hobbs et al. 1985] Hobbs L.M., Vidal-Madjar A., Ferlet R., Albert C.E. & Gry 
C., 1985, ApJ, 293, L 29 

[Holland et al. 2003] Holland W.S., Greaves J.S., Dent W.R.F, et al., 2003, ApJ, 
582, 1141 



12 



[Jayawardhana et al. 2002] Jayawardhana R., Holland W.S., Kalas P., et al., 2002, 
ApJ, 570, L 93 

[Kalas & Jewitt 1995] Kalas P. & Jewitt D., 1995, AJ, 1 10, 794 

[Kamp et al. 2003] Kamp I., van Zadelhoff G.J., van Dishoeck E.F. & Stark R., 
2003, A&A, 397, 1129 

[Kuchner & Holman 2003] Kuchner M.J. & Holman M.J., 2003, ApJ, 588, 1 1 10 

[Lagrange et al. 2000] Lagrange A.-M., Backman D.E. & Artymowicz P., 2000, 
in: Mannings V., Boss A.P Russell S.S. (eds.) Protostars and Planets IV, Uni- 
versity of Arizona, p. 639 

[Lecavelier des Etangs et al. 2001] Lecavelier des Etangs A., Vidal-Madjar A., 
Roberge A., et al., 2001, Nature, 412, 706 

[Liseau 1999] Liseau R., 1999, A&A, 348, 133 

[Olofsson et al. 2001] Olofsson G., Liseau R. & Brandeker A., 2001, ApJ, 563, 
L77 

[Ozernoy et al. 2000] Ozernoy L.M., Gorkavyi N.N., Mather J.C. & Taidakova 
T.A., 2000, ApJ, 537, L 147 

[Pantin et al. 2000] Pantin E., Els S., Marchis E, et al., 2000, AAS-DPS 32, 
#65.33 

[Roberge et al. 2000] Roberge A., Feldman P.D., Lagrange A.M., et al., 2000, 
ApJ, 538, 904 

[Schneider et al. 1999] Schneider G., Smith B.A., Becklin E.E., et al., 1999, ApJ, 
513, L 127 

[Schneider et al. 2001] Schneider G., Becklin E.E., Smith B.A., et al., 2001, AJ, 
121,525 

[Smith & Terrile 1984] Smith BA. & Terrile R.J., 1984, Science, 226, 1421 

[Spangler et al. 2001] Spangler C., Sargent A.I., Silverstone M.D., Becklin E.E., 
& Zuckerman B., 2001, ApJ, 555, 932 



13 



[Takeuchi & Artymowicz 2001] Takeuchi T. & Artymowicz P., 2001, ApJ, 557, 
990 

[Thiet al. 2001] Thi W.F., Blake G.A., van Dishoeck E.F., et al., 2001, Nature, 
409, 60 

[Trilling et al. 2000] Trilling D.E., Brown R.H. & Rivkin A.S., 2000, ApJ, 529, 
499 

[Wahhaj et al. 2003] Wahhaj Z., Koerner D.W., Ressler M.E., et al., 2003, ApJ, 
584, L 27 

[Wilner et al. 2002] Wilner D.J., Holman M.J., Kuchner M.J. & Ho P.T.P., 2002, 
ApJ, 569, L 115 

[Wyatt & Dent 2002] Wyatt M.C. & Dent W.R.F., 2002, MNRAS, 334, 589 
[Zuckerman 2001] Zuckerman B., 2001, ARAA, 39, 549 



14 



Figure Captions 



Fig. 1. The ratio of integrated excess to photo spheric flux as a function of 
estimates of the stellar age. The most prominent Vega-excess stars are identified. 
Open circles refer to stellar members of a number of open clusters. The point for 
the Sun is a lower limit, as it traces only the fraction of the zodiacal dust inside 
the Kuiper Belt. Adapted from Spangler et al. (2001). 

Fig. 2. Observations of the (3 Pic disk in the Nal D 2 A5890 line (ESO, Chile). 
In the right panel of the Left Figure, the EMMI {R ~ 60 000) observations with 
the 3.5 m NTT are shown (horizontal scale in kms -1 , vertical scale in AU). A 
model of the line emission from a disk in Keplerian rotation is shown to the left, 
whereas, in the middle, the model is shown degraded to the resolution of the 
observations, with w hite noise added and the central stellar spectrum removed 
( Olofsson etal. 200 ll) . Follow-up observations with UVES (R ~ 90 000) at the 
8 m VLT, shown in the Right Figure to scale, largely improve d on both the spec- 
tral and spatial resolution and coverage dBrandeker et al. 20031) . Telluric D 2 emis- 
sion is seen at — 30kms -1 and the approaching southwest side is up (cf. Fig. 3). 

Fig. 3. The slit positions of our VLT-UVES observations are shown super- 
posed onto the HST-STIS image of the scattering dust disk ( Heap et al. 2000h 
projected onto the sky. The width of the slit is 0"3 and, during the observations, 
the seeing was in the range 0'.'4 to 0'.'6. The length of the slit, limited by the Echelle 
inter-order separation, is 8". The observations along the disk were performed with 
overlapping slit positions. 

Fig. 4. Scale heights of the gas disk in three emission lines (ground state tran- 
sitions). The bars show the Full Width to Half Maximum of the line emission 
at four positions in the disk, x- and ^/-offsets, in arcs ec, are relative to the star 
j3 Pic and to the midplane of the scatterin g dust disk dKalas & Jewitt 19951) . re- 
spectively. For comparison, the dust disk (IHeap et al. 200 0) is shown to scale in 
the lower half of the figure. 

Fig. 5. The spatial distribution of the NalD? emission in the j3 Pic disk, with 
the NE and SW disk halves identified dBrandeker et al. 20031) . Error bars represent 
la of the quadratically summed photon noise of the line and off-line fluxes. Radial 
distances from the star are expressed in arcsec and integrated surface intensities 
in erg cm~ 2 s _1 arcsec -2 . The vertical dashed lines refer to the positions of the 
perpendicular slits. 

Fig. 6. The UV spectrum of (3 Pic. Continuum points are from FUSE and HST 
observations (archive data). Also shown are the results of HST-STIS observations 
of the NE- and SW-disk of j3 Pic at 0"6 distance from the star (archive data). 



15 



The line spectrum is based on a stellar m odel atmosphere (at 19 .3 pc) with T cS = 
8000 K, log g = 4.5 and log Z/Z & = 0.0 dHauschildt et al. 1999h : for the assumed 
R = 1.75 R Q , the luminosity of the star is L = 11 L Q . The ionisation edges of 
species discussed in the text are indicated by the vertical bars. 

Fig. 7. The asymptotic velocity versus the radial distance r from the star. 
For Nal (A = 23) an average [3 = 250 has been used, and the red curves are for 
n(26 AU) = 10 6 , 10 7 and 10 8 cm^ 3 forp = 0, -1.5 and - 2.5 (dotted, dashed 
and solid lines) respectively (see the text). For Call (A = 40) /3 m j n = 3 has been 
applied as a strict lower limit to the radiation pressure coefficient (see Table[T]). 
The blue curves are for n(26 AU) = 10 4 ' 5 cm^ 3 for p = and 10 5 cm -3 for both 
p — —1.5 and —2.5. For QTI (A = 52), only the constant density case (10 5 cm" 3 ) 
is shown by the green curve. As this ion is insensitive to radiation pressure, rela- 
tive velocities are negative. The 5a upper limit to the observed excess velocities 
over Keplerian motion is shown by the purple horizontal line. 
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